Uniform Semiclassical Approach to Fidelity Decay in the Deep Lyapunov Regime 
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We use the uniform semiclassical approximation in order to derive the fidelity decay in the regime 
of large perturbations. Numerical computations are presented which agree with our theoretical 
predictions. Moreover, our theory allows us to explain previous findings, such as the deviation from 
the Lyapunov decay rate in cases where the classical finite-time instability is nonuniform in phase 
space. 
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The stability of quantum motion under a system's per- 
turbation can be measured by the so-called quantum 
Loschmidt echo or fidelity Q, 0, El • It is defined as the 
overlap M(t) — \m(t)\ 2 of two states obtained by evolv- 
ing the same initial state |$o) under two slightly different 
Hamiltonians: 



m{t) = (* |exp(ii«/ ? i)exp(-iffot/ft)|*o) 



(1) 



Here Hq is the Hamiltonian of a classically chaotic sys- 
tem and H — Hq + eV is the perturbed Hamiltonian, 
with e a small quantity and V a generic perturbing po- 
tential. This quantity can also be seen as a measure of 
the accuracy to which an initial quantum state can be 
recovered by inverting, at time t, the dynamics with the 
perturbed Hamiltonian H . 

This quantity has attracted much attention recently, 
mainly in relation to the field of quantum computation 
and in connection to the corresponding classical motion 
HSEIHBUHEIEIE!- Focusing on systems with 
chaotic classical limit, one may identify, by increasing the 
perturbation strength, three different regimes of fidelity 
decay: (i) The perturbative regime, in which the fidelity 
has a Gaussian decay, (ii) The Fermi-golden-rule regime, 
with an exponential decay of fidelity, M(t) oc exp(— Ft). 
Here the decay rate T is the half-width of the local spec- 
tral density of states (LDOS) 0, which can also be cal- 
culated semiclassically 0. (iii) The Lyapunov regime, in 
which M(t) oc exp(— At), with A being the (maximum) 
Lyapunov exponent of the underlying classical dynamics 

a- 

However, the above picture remains unsatisfactory. 
This is particularly the case in the deep Lyapunov regime 
with a 3> 1, where a — e/h is the parameter character- 
izing the strength of quantum perturbation, and Ti the 
(effective) Planck constant. In systems with nonconstant 
finite-time Lyapunov exponent (which is the typical sit- 
uation), fidelity decays with a rate different from A [llj . 



Indeed, a semiclassical analysis ^l] leads to an exponen- 
tial decay of fidelity with a rate Ai < A. The relation 
between this semiclassical treatment and that along the 
lines of Refs. 0,|EGlH>Ei is unclear. Moreover, an ex- 
tremely fast, super-exponential decay of fidelity has been 
found within a quite short initial time for initial Gaussian 
wave packets |l 1| . In view of the importance of fidelity for 
the characterization of the stability of quantum motion 
under a system's perturbation, it is necessary to provide 
a clear theoretical understanding of its behavior and, in 
particular, to account for the seemingly disconnected and 
sometimes contradictory results. 

In this paper, we focus on the behavior of fidelity 
for a > 1 and we treat this problem in full general- 
ity. We derive the general semiclassical formula which 
correctly reproduces the two limiting cases of exp(— At) 
and exp(— Ait) decays. We also show that under certain 
conditions the exponential rate of fidelity decay can be 
equal to twice the classical Lyapunov exponent. 

Our starting point is the semiclassical approximation 
to the fidelity for an initial Gaussian wave packet given 
in El, 



m sc (t) ~ (£ 2 M ) 



2\<V 2 



J rfpo exp [iA5( Po , r ; t)/h - (p - po) a /(h/t) 3 



(2) 



where AS'(po, ro! t) is the action difference along the two 
nearby trajectories starting at (po, ro) m the two systems 
H and Ho, 

AS( Po ,?o;t)-e [ dt'V[r'(t')] (3) 



with V evaluated along the trajectory in the Hq system. 
The initial Gaussian wave packet, centered at (ro,po)j is 



^o(ro) = K 2 ) ^ exp [ ? p • r /h - (r - ? ) 2 /(2£ 2 



(4) 
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For simplicity, we will consider here kicked systems with 
d = 1 and set the domains of r and p to be [0, 2ir). The 
effective Planck constant is taken as h — 2ir/N, where N 
is the dimension of the Hilbert space. 

The main feature of AS as a function of po is its oscil- 
lations, the number of which increases exponentially with 
time t. Indeed, the variance of AS increases linearly with 
t 0, while the slope of AS/e, denoted by k p , 



_ 1 dAS( Po ,r ;t) 
P e dp a 



* ,dV dr'(t') 
o dr' dp 



(5) 



increases on average exponentially with t, due to the fac- 
tor dr' /dpo. 

Let us first discuss the fidelity for a single initial state. 
Neglecting a quite short initial time, the main contri- 
bution to the right hand side of Eq. (J2J comes from 
the integration over the region [po — w p ,po + w p ], where 
w p = h/£ is the width of the initial Gaussian in the po 
space. Let us define a time scale r such that at t = r, 
AS(pQ,r ;t) completes one full oscillation period as po 
runs over [po — w p ,po + w p ]. In a system possessing a 
constant local Lyapunov exponent A, the number of os- 
cillations of AS increases exponentially as coe At and we 
obtain the average estimate for t, 



1 



\h(tt/c w p ). 



(0) 



Taking, e.g., £ = \f%, one can clearly see that r is of the 
order of the Ehrenfest time j \nfi~ 1 . In the general case 
of systems with fluctuation in the finite-time Lyapunov 
exponent, the number of the oscillations of AS increases 
as e A (*)*, with some time-dependent rate A(i) < A. 

We consider first the behavior of fidelity for t < t and 
denote with k p the value of k p in the center (po,ro) of 
the initial Gaussian. For such times, the phase AS/h, on 
the right hand side of Eq. @, as a function of po, can 
usually be approximated by a straight line with a slope 
ak p , within the region po € [po - w p ,po+w p ]. Due to both 
the fast increasing of \k p \ with time and the large er value, 
one has \ak p \ 3> n/w p for most initial states and, as a 
result, the change of the phase (AS/h) within the interval 
Po € [po — w p , po + w p ] is much larger than 2ir. Note that 
the largest slope of the term (po — Po) 2 /w p within this 

interval of po is 2/w p , which is much smaller than |<rfc p |. 
The right hand side of Eq. @ can now be calculated 
approximately within the interval po S [po — w pi po + w p ] 
and gives 



M SG (i) tx l/{akp) 2 



(7) 



For times t > r, or when \k p \ is small enough for t < r, 
the stationary phase approximation can be used in calcu- 
lating m sc (t) in Eq. If we denote by a the stationary 
points and by po a the momenta at which k p = 0, we have 
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FIG. 1: Decay of averaged fidelity in the map 1151 . with 
K = 1, rj = 0.987. I s and Ia are the theoretical predictions 
1 1 21 and 1141 1. respectively. It is seen that after a short initial 
time, both I s (t) and I\(t) are close to the exact fidelity M(t) 
(until saturation is reached). For comparison, the decay e _Alt 
is shown. We also plot the average classical fidelity, which was 
calculated by taking initial points within circles with radius 
Vfo in the phase space [Tct fl^| . For this map, A w 0.9 and 
Ai w 0.35. Here and in the following figures, a — 100, N — 
131 072, £ = Vh, and averages are performed over 2000 initial 
Gaussian packets. 



J2 a m a (t), where 
2ihexp [j-AS(p 0a ,¥ ;t) - (p 0a - p ) 2 /Wp] 




with AS„ = 



d 2 AS( Po ,r ;t) 



dpi 



(8) 



Next we turn to the behavior of average fidelity and 
first consider the long time decay, namely, t > t. Due to 
the large a value and to the classically chaotic motion, 
the phase AS(po a ,ro;t)/h in Eq. © can be regarded as 
random with respect to a and 7q. Then, the averaged 
fidelity M(t), with average taken over Fo and po, can be 
approximated by its diagonal part 01 > 



(9) 



The right hand side of Eq. © can be expressed as an 
integration of l/|fc p |. For this, we introduce A a to denote 
the region [Po a ,Po a -S] U[Poa+£,PoJ, where pg Q = {p 0a + 
po,a-i)/2, Poa = (poa +po,a+i)/2, and where 6 is a small 
quantity. In the neighborhood of po Q , k p satisfies ek p ~ 
ASq, • (po — poa)- For small enough S, we have 



dp - 



ck v 



2 In 5 



(10) 



Substituting the expression of |AS^| obtained from 
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FIG. 2: Same as Fig.0 but for 77 = 0.85, for which A « 0.92 
and Ai « 0.81. 



Eq. (TJJ i nto Eq. ©, we have 



M(t) 



(27r) 3 /2(_21n5)e 



dr / dp - 



A',, 



(11) 



where := LL^*- 

Since the value of 6 is irrelevant for the decay rate, we 
may write 



M(t)ocI s (t) := 



'/'■<) / ^0 77^- 



(12) 



An accurate numerical evaluation of I s (t) is not easy 
since one must find out all stationary points a for each 
value of r . An approximate numerical result can be 
obtained by using the Monte Carlo method in which, in 
order to perform the integral (|12f> over the region Vg, i.e., 
with the neighborhoods of stationary points excluded, we 
neglect the small set of points that have the smallest val- 
ues of \k p \. 

Actually, one can make a further approximation by 
using the following arguments. The main contribution 
to the integral in Eq. I|ll|) comes from small values of 
\k p \ in the region Vg. For p £ Vg close to a stationary 
point p Qa , k p in Eq. (J5J can be approximated by 



dt' 



d 2 v / a 



dr 12 \ dpo 



dr' dp 2 , 



(P0-P0a)- (13) 



Due to exponential divergence of neighboring trajecto- 
ries in phase space, the main contribution to the right 
hand side of Eq. fTSjl comes from times t' as t. The time 
evolution of the quantity inside the bracket in Eq. i|13|) 
is given by the dynamics of the system described by Hq. 
On average it increases as [Sx(t)/5x(0)] 2 , where 5(x) de- 
notes distance in phase space. With increasing time, the 
number of the stationary points of AS increases exponen- 
tially, roughly in the same way as Sx(t)/5x(0), since the 
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FIG. 3: Decay of averaged fidelity in the sawtooth map (77 = 
0) with K = 1 and i = 3 in Eqs. and ifT?). showing e~ 2Xt 
decay followed by the Lyapunov decay (A = 0.96). The arrow 
indicates the theoretical estimate of the crossover time r. 



oscillation of AS is mainly induced by local instability of 
trajectories. Then, substituting Eq. i|13|) into Eq. I|ll|l . 
we have M(t) oc |<5a;(t)/fe(0)| -1 , which can be written 
as 

M(t) oc I A (t) = exp[-Ai(*)i], 
1 



with 



Ai(i) 



lim In 

t <5a;(0)^0 



Sx(t) 



5x(0) 



(14) 



In systems with constant local Lyapunov exponents, 
Eq. (|14fl reduces to the usual Lyapunov decay with 
Ai(t) = A. On the other hand, when fluctuations in 
local Lyapunov exponent cannot be neglected, I/v(i) co- 
incides with the e~ Xlt decay in Ref. 01 with Ai = 
limt_>oo Ax (t) < A, only in the limit t — > 00. Therefore, 
the actual decay, which can be observed in finite times, 
can be considerably different from the e~ Xlt decay. 

For times t < r, the main contribution to the aver- 
aged fidelity M(t) comes, in fact, from initial states with 
k p close to zero. When AS is not too flat, one can still 
use the stationary phase approximation for these initial 
states. Hence, also for i<rwe obtain the same expres- 
sions as in Eqs. (|12fl and l|14|) for the decay of averaged 
fidelity. 

We have tested the above predictions by considering 
the map 

Pn+i = Pn + K[{r n -7r) +7jsinr n ] (mod 2n), 
r n +i = r n +p n+1 (mod 2?r), (15) 

with two parameters K,rj £ [0, 1]. For K > and r\ = 0, 
this is the piece- wise linear sawtooth map |lfjj . which is 
hyperbolic with constant local (finite time) Lyapunov ex- 
ponent. For the particular case K — 1, the map reduces 
to the perturbed cat map, which is known to be Anosov 
for < r] < 1 (having non-constant A), whereas for rj = 1 
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FIG. 4: Fidelity decay in the sawtooth map (7/ = 0) with 
K — 1 and i — 2. Left panel: M(t) of a single initial Gaussian, 
showing large fluctuation at t < 4, approximate e~ 2A * decay 
within 4 < r: < 7, and approximate Lyapunov decay at t > 
8 (before saturation). Right panel: averaged fidelity M(t), 
showing the Lyapunov decay. 

it acquires a marginally stable (parabolic) fixed point. 
This map is quantized in a Hilbert space of dimension N. 
The one period quantum evolution is given by the Flo- 
quet operator, U = exp[— ip 2 /(2fr)] exp[— iU(r)/%], with 
U(r) = — K[(r — ?r) 2 /2 — 7ycosr]. In order to compute 
fidelity, we choose to perturb the parameter K — > K + e. 
Figure ^ shows that numerical data accurately fit our 
theoretical predictions in Eqs. (fL2"|) and l[Tl)l. In Fig. 
it is seen that with decreasing 77, since the values of A 
and Ai become closer, our predictions approach that of 
Ref. [llj. At 77 = 0, the classical map has a constant local 
Lyapunov exponent and the standard Lyapunov decay is 
recovered. 

In the above discussion of the average fidelity, the exis- 
tence of stationary phase is assumed. It may happen, in 
some circumstance, e.g., with some special perturbation, 
that there is no stationary phase for AS. In this case it 
turns out that a decay with a rate of double Lyapunov 
exponent may appear for t < r, when the classical sys- 
tem has a constant local Lyapunov exponent. Indeed, 
for t < t, the main contribution to the averaged fidelity 
M{t) comes from initial states for which the values of 
\k p \ are close to local minimum of \k p \. When the values 
of local minimum of \k p \ are large enough, the decay of 
the fidelity is given by Eq. Then, since \k p \ increases 
on average as e At , the averaged fidelity has a double- 
Lyapunov-exponent decay, 

M(t) oc e~ 2A \ t < t. (16) 

For t > r, one can use arguments given in Ref. fl3| . 
showing that M(t) still follows the standard Lyapunov 
decay, M(t) oc e~ At . 

Finally, in systems possessing stationary phase in AS" 
and constant local Lyapunov exponents, although the av- 



eraged fidelity has Lyapunov decay, a double-Lyapunov- 
exponent decay e~ 2At may appear for t < r, for the fi- 
delity of those single initial states, for which \ k p \ happens 
to increase exponentially as e At [see Eq. Q]. 

In order to check the above predictions, we consider 
the sawtooth map (7/ = 0) which has a constant local 
Lyapunov exponent, A = ln({2+K +[(2+K) 2 -4} 1 / 2 } /2). 
We consider here the following perturbed map 

Pn+i = Pn + K(r n - 7r) + eij\fi(r„ - 7r) l_1 , i = 2, 3, 
r n +i = r n +p n +i, (17) 

with A/2 = 1/2 and A3 = VTA/3tt. These two values A/j 
give the same decay rate in the Fermi-golden-rule regime. 
However, while for i = 2 stationary phase of AS exists, 
in the case i = 3 there is no stationary phase in AS* vs po- 
In the latter case, as shown in Fig. [31 the average fidelity 
has an initial double-Lyapunov-exponent decay followed 
by the standard Lyapunov decay, as predicted by the 
theory. The crossover of the two decays is in agreement 
with the theoretical estimate r ~ 6.5. Figure 01 (left 
panel) shows instead that a double-Lyapunov-exponent 
decay may appear for the fidelity of some particular single 
initial state, while the average fidelity has the Lyapunov 
decay (right panel). 

In summary, we have derived general semiclassical ex- 
pressions for the fidelity decay, at strong perturbations, 
which reproduce, as two particular limiting cases, pre- 
vious results leading to the Lyapunov decay and to the 
e -Ait decay. In particular we have discussed the rele- 
vance of fluctuations in the finite-time Lyapunov expo- 
nent and we have shown that fidelity decay depends on 
the strength of such fluctuations in the Lyapunov regime. 
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